Does Dirac-Born-Infeld Modification Of Quadratic Theories Really Matter? 
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We study the consequences of further modification of f{R, R^^R'^" , Rf^vapR'^""'') / f{R)-t]ieones 
by means of th e Dirac-Born-Infeld deformation procedure, which amounts to the replacement of / by 
X{^/T^\-'2JJX — 1) (the free parameter A fixes an additional energy scale). We pay special attention 
to the definition of masses of the linearized propagating degrees of freedom since these are important 
to judge about the stability of the linearization around vacuum background spaces. In this context 
we discuss the subtleties associated with expanding f{R, Rui^R'^" , Rfj,va pR'^"" '')-La,gTa,ngia,ns around 
maximally symmetric spaces of constant curvature, as well as with equivalence of the linearized 
Lagrangian to a scalar-tensor theory. Investigation of the consequences of applying the Dirac-Born- 
Infeld strategy to further modify quadratic theories, on the stability of de Sitter vacuum, as well as 
its impact on the cosmological dynamics, is the main concern of this paper. 

PACS numbers: 04.20.-q, 04.50.Kd, 95.36.-l-x, 98.80.-k, 98. 80. Bp, 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 



Attempts to modify the Einstein-Hilbert (EH) action 
of general relativity (GR) 

Seh^^J d^x^\{R~ 2A), 

where R = g'^'" R^v is the Ricci curvature scalar, and A- 
the cosmological constant (k^ = rUpf = 8nG) , have been 
motivated by a number of reasons. In particular, renor- 
malization at one-loop demands that the Einstein-Hilbert 
action be supplemented by higher order curvature terms 
Besides, when quantum corrections or string theory 
are taken into account, the effective low energy action for 
pure gravity admits higher order curvature invariants Q . 

More recently it has been suggested that the present 
cosmic speed-up could have its origin in - among other 
possibilities - corrections to the GR equations of motion, 
generated by non-linear contributions of the scalar cur- 
vature R in the pure gravity Lagrangian of f{R) theories 
0-0] ■ Solar system constraints on f{R) theories that 
are able to accommodate present accelerated expansion 
of the Universe, have been one of the most discussed 
subjects lately |8l-[l3|. see also [13] for an extensive re- 
view. Comparison of these theories with solar system 
measurements, relies on the weak-field limit expansion of 
the f{R) Lagrangian, and the consequent calculation of 
post-Newtonian contributions to the metric coefficients 
[1, Slill- Nonetheless, even if f{R) theories were not a 
viable alternative to explain current acceleration of the 
expansion, their relevance to study early-time infiation 
[15| might fuel further interest in these alternatives to 
general relativity. 



Next in degree of complexity are the so called 
f{R,Rfjj^>'\Rf,^apR'''''''') - fiR,-) for short - theo- 
ries [3, HB-Hol ■ The gravitational spectrum of the lin- 
earization of these theories consists of a massless spin-2 
graviton plus two additional gravitational propagating 
degrees of freedom: a massive spin-0 excitation, and a 
massive spin-2 propagating mode. The latter appears 
to be a ghost mode associated with the Weyl curvature 



invariant = Cn 



ucrX 



Notwithstanding, 



there are ways to overcome (or at least to smooth out) 
the consequences of the would be massive spin-2 ghost 
mode [ill. 

There are additional ways to modify the EH GR ac- 
tion. For instance, the one based on the Dirac-Born- 
Infeld (DBI) procedure for smoothing out singularities 
[3 mi [13 According to this procedure the original La- 
grangian density C = y/\g\L - whose singularities are to 
be cured - is replaced by one of the DBI form: 



C 



C 



DBI 




where the scale A sets an upper bound to curvatures ac- 
cessible to the theory. A combination of the above possi- 
ble modifications, i. e., a DBI- type action containing an 
f{R, ...) function within the square root, i. e., a replace- 
ment of 



/(E,...)^ A 




2/(i?,.. 
A 



(1) 



in the action S cc J d'^x^y\g\f{R,...), could supply an 
additional cosmological scenario where to look for altcr- 
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^ Higher order actions are indeed renormalizable (but not unitary) 
i). 



^ The proposal to remove initial as well as final singularities in 
modified gravity has been given in Ref. l23l . It was shown there 
that the addition of a term to otherwise divergent modified 
gravity makes it regular. 
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native explanations to several phenomena such as infla- 
tion and the present speedup of the cosmic expansion. 
Several theories of gravity of this kind have been pro- 



posed since long ago in [24| . and in more recent years, 
for instance, in [2^ (see also [1^). To be phenomenolog- 
ically viable, non-linear modifications of general relativ- 
ity have to satisfy several physically motivated require- 
ments [1^ : i) reduction to EH action at small curvature, 
ii) ghost freedom, iii) regularization of some singularities 
(as, for instance, the Coulomb-like Schwarzschild singu- 
larity), and iv) supersymmetrizability. Nonetheless, the 
latter requirement is quite stringent and, for most pur- 
poses, might be excluded. 

Would further modification of f{R, ....)/ f{R) theories 
of gravity, through the DBI deformation strategy ([T]), 
help surmounting the severe problems related with the 
presence of a multitude of instabilities within these the- 
ories? Would it modify the asymptotic properties of the 
cosmic dynamics? Aim of the present paper is, precisely, 
to investigate the consequences of applying the Dirac- 
Born-Infeld procedure to /(i?, ...) theories of gravity - 
including /(i?) gravity as a particular case -, regarding 
stability of de Sitter vacuum solutions, as well as its im- 
pact on the cosmological dynamics. 

The paper has been organized as follows. In section II 
the subtleties associated with expanding f{R, ...) gravity 
theories about maximally symmetric spaces of constant 
curvature are discussed. This topic is central to judge 
about stability of the propagating degrees of freedom 
upon linearization. The equivalence of ghost-free f{R, ...) 
gravity to a scalar-tensor theory is demonstrated in the 
same section. Sections III and IV are devoted to study 
modifications of the stability properties of f{R, ...), and 
f{R) theories, respectively, after applying to them the 
DBI procedure. The consequences for the cosmological 
dynamics of DBI-modified f{R,...)/f{R) gravity is the 
main concern of section V. The results of the present 
investigation are discussed in section VI, while the con- 
clusions are given in section VII. An appendix with brief 
and concise tips on how to apply the dynamical systems 
tools, is included in the final section VIII. 



II. EXPANSION AROUND MAXIMALLY 
SYMMETRIC SPACES OF CONSTANT 
CURVATURE 



where R is the curvature scalar, R^^, - the Ricci ten- 
sor, and Rp,ua\ - the Riemann tensor. In terms of these 
variables the Gauss-Bonnet invariant Q and the Weyl in- 
variant = C^jyi^AC'^'^'^^ can be written as 



g = - AY + Z , and = Z - 2r + -X^ , (3) 

3 

respectively. Combining these expressions one can get 
the following equalities that will be useful latter on; 



Y = -c'^ --g + -x^ . z = 2c^ -g + -x^ . u) 

2 2 3' 3 ^ ' 

We will consider pure gravitational actions of the fol- 
lowing kind: 



Sn 



1 

2^ 



dxW\g\ f{X,Y,Z) 



(5) 



where = SttGn = Mpf. The following field equations 
can be obtained from the above action by varying with 
respect to the metric g^i, 0: 



fxG^, = -g^^if - Xfx) - [g^un - V^V,)/x 

-2{fYR%Ra. + fzRxap^R^^",) - 9^.^ xifv R"^) 

~U{fyR,,) + 2V.Va (/yi?^(^i?\) + 2/zi?^(^,)^) ,(6) 

where, as usual, G^,y = — g^^R/2 is the Einstein 
tensor, □ = g^'^V ^V^, is the D'Lambertian, and we 
have used the usual representation for symmetrization: 
T'(tiv) = -\- T^fj)/2. In the presence of matter these 
field equations amount to: 



G'^.-8^G^/^'(r^r' + ?;T), 



(7) 



where ^ttG"^^ - 



K / fx is the effective gravitational cou- 
pling, T^™'' is the stress-energy tensor for matter, while 
T™'" equals times the right-hand-side (RHS) of equa- 
tion ([6]). The trace of equation (|6]) generates an addi- 
tional constraint on the curvature: 



Here we try exposing the importance of considering a 
consistent expansion around maximally symmetric vac- 
uum spaces of constant curvature of higher order f{R, ...) 
theories, to judge about stability issues. The results ob- 
tained can be applied also to f{R) theory as a particular 
case. In what follows, for simplicity of writing we use the 
following definition of variables X^ — (X, y, Z): 



X = R , Y = RnuR^'^ , Z = Rnv(y\R^ 



(2) 



2/ - Xfx - 2YfY - 2Zfz 
-□(3/x + 2Xfy + Xfz + Ix^fy) = A^'T(") (8) 

Let us consider expanding the above action ^ around 
maximally symmetric vacuum spaces of constant curva- 
ture R — Rq, i. e., we will Taylor expand f{X,Y,Z) in 
the neighborhood of the point {Xo,Yo, Zq) where Xq = 
Ro, Yo = Rl/A and Zo = Rl/6, and f{Xo,Yo,Zo) = /o, 
up to the second order ~ {X ~ Xq)'^ ^ (Y — Yq)^ ^ 
{Z - Zof. One has: 



3 



fiX,Y,Z) = ,fo + Y,(-^) {X"-X^) 



1 



2 \dX'dX^ 
or, in explicit form: 



{X^ - Xl,){X'' - X^) + 0{i) ,(9) 



f{X, r, Z)^h + fx{X - Xo) + f^{Y - Yo) 

+f°{Z-Zo) + lfxxiX-Xor 
+f^xiX - Xo){Y - Yo) + fzxiX - Xo){Z - Za) 
+ lfYYiY - Yof + fMY - YoKZ - Zo) 

+ \fzz{Z-Zof+0{i), (10) 



where 



Jx 



dX 



J XX 



\dXdxJ,' 



etc. Notice that, while keeping up to the second order in 
the expansion is legitimate since we are considering small 
SX"^ — X'^ — Xq <C 1, considering up to quadratic terms in 
the curvature ^ is an additional requirement that has 
nothing to do with the order of the perturbations around 
Xq = {Xq, Yo,Zo) one is being considering. Actually, the 
expansion (I9ll0p does not exclude large curvature Rq. 
The only requirement is that if Xq — Rq is large, so is X. 
However, in this latter case it makes sense to keep only 
the terms with the higher orders in the curvature (terms 
ocX\Y\Z^). 

After had clarified this point let us consider small cur- 
vature backgrounds only, i. e., small Rq, so that it makes 
sense to keep terms up to the order ~ 0{R^) in (fTO)) : 



As a check of consistency, notice that by taking into 
account the trace equation ([5]), which in the present sit- 
uation amounts to: 



2/o - Xofx - 2Yof^ - 2Zofz = , (14) 
then the coefficient 

^0 = — 2^*--^^ ^ ^ofxx) = — ^-'^o ■ 

This means that the definition of the cosmological con- 
stant A = — Ao/2ao (see below) leads to the expected 
result A = Ao/4. 

As customary one can introduce, additionally, the fol- 
lowing magnitudes: 



ao 



Jx ^ ^oJx 



XX 



70 



-2A=^ 



so that (IT2]) can be written in the following form: 



(15) 



/(X,y,Z) = ao -2A + i? 



1 i?2_ 1 ^2 



6mQ 



2ml 



(16) 

The gravitational spectrum of the corresponding lin- 
earized theory consists of a standard massless spin-2 ex- 
citation, plus a spin-0 mode of mass squared mi , and 
an additional spin-2 mode of mass squared m\ [la, U^l - 
TOq and ml are given by p^ . The latter spin-2 excita- 
tion occurs to be a ghost, leading to non- unitary states 
upon quantization (otherwise unitary states of negative 
energy). ^ 



f{X, Y, Z) — ff) - Xofx - Yofy ^ ^ofz 

+ 2^ofxx + {fx - ^nfxx)^ 

+ \fxxX^ + f^Y + fzZ + 0{R^) .(11) 

By using the equations Q one can remove the vari- 
ables Y , Z from the expansion, so that, finally: 



/(X,y,Z) = Ao + ao^+^A2 + ^C2 , (12) 



where 



A. Stability Requirements 

Stability issues are central in the study of higher or- 
der modifications of general relativity, since these are 
plagued by several kinds of instabilities, some of which 
are catastrophic, leading to subsequent ruling out of the 
corresponding theories. Amongst these is the fundamen- 
tal Ostrogradski instability, based on the powerful no-go 
theorem of the same name 27[ : "there is a linear insta- 
bility in the Hamiltonians associated with Lagrangians 
which depend upon more than one time derivative in 
such a way that the dependence cannot be eliminated 
by partial integration" . This result is general and can 
be extended to higher order derivatives in general. As 



X^ 

Ao = /o - Xo!x ^ Yofy - Zofz + —fxx > 

"0 = /x - Xofxx 1 = ^fxx + 2/y + 2fz 

70^/0+4/0. (13) 



^ In Ref. [Tsll this has been shown starting directly from a general 
action of the form (|5)l by investigating the propagator of the 
linearized degrees of freedom. 
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a consequence, the only Ostrogradski-stable higher order 
modifications of Einstein-Hilbert action are those in the 
form of an f{X) = f{R) function [l^l- This result alone 
might rule out any intent to consider higher order mod- 
ifications such as quadratic ones, for instance. However, 
the subject is subtle and, in the last instance, consider- 
ation of such theories can shed more light on the stabil- 
ity issue. As an example, consider theories of the kind 
f{X,Y,Z) where the invariants Y, Z, enter in the same 
combination as in the Gauss-Bonnet invariant 5, i. e., in 
the combination —4Y + Z. In this case, upon linearizing 
around maximally symmetric spaces of constant curva- 
ture, since fy +4/^ = 0, the massive propagating spin-2 
ghost associated with the Weyl invariant decouples from 
the gravitational spectrum. Besides, the linearized the- 
ory can be recast into the form of an /(X)-theory (see 
Eq. (IT6|) '). so that it is also Ostrogradski-stable. 

Another kind of catastrophic instability is the so called 
"Ricci instability" [JJj], also known as Dolgov-Kawasaki 
instability [28]. For arbitrary /(X)-theories the analysis 
of Ricci stalaility has been generalized in [2^ , while a con- 
sistent and simple physical interpretation has been given 
by the same author in Ref. [s^l- According to the latter 
interpretation, assuming that the effective gravitational 
coupling SwCj^-^ = k'^/ fx (see Eq. ([7])) is positive, then, 
since 



dG 



eff 

N 



dX 



for negative fxx < 0, the effective gravitational cou- 
pling G^'^ increases with the curvature, otherwise, at 
large curvature gravity becomes stronger, and since X 
itself generates larger and larger curvature through the 
trace equation Eq. ([8|), the effect becomes uncontrol- 
lably stronger because of an increased C)^^ . In other 
words, a positive feedback mechanism acts to destabilize 
the theory [llj. Hence, to avoid the "Ricci instability" it 
is necessary that fxx > 0. 

Other instabilities, such as those caused by the pres- 
ence of a spin-0 tachyon degree of freedom are not of less 
importance. The latter instability is associated with neg- 
ative values of the spin-0 mass squared in Eq. ([T5|) . 
Last but not least, there is an additional requirement 
that has not been discussed in detail in the literature. 
As seen from p6)) . the constant = f^ — Xof^-^ is 
an overall factor that multiplies the linearized action ([5]) , 
and, hence, it may change the sign of the action, other- 
wise, the sign of the effective gravitational coupling upon 
linearization. Therefore, if one considers f{X) theories; 
Ricci stability, absence of tachyon, and positivity of the 
effective gravitational coupling, are not independent re- 
quirements. Actually, in this case 



nin 



fO 

Jx 



^ofxx 



XX 



so that, Ricci stability {fxx > 0), and positivity of 
the effective gravitational coupling upon linearization 



(fx — Xofxx > 0)j together imply absence of tachyon 
instability ttzq > 0. For arbitrary f{X, Y, Z) the above 
requirements are independent. 

To summarize the discussion on the relevant stability 
requirements imposed on higher order curvature modifi- 
cations of Einstein-Hilbert theory, here we list them: 

• Ostrogradski Stability: The linearized f{X,Y,Z) 
theory should be expressible as an equivalent f{X)- 
theory.^ This can be implemented if 



/o +4/1-0 



• Ricci Stability: 



fxx > 



• Absence of Tachyon Instability: 

™2 _ Jx ^ ^oJxx 



mx+Vy + Vz 



> 



• Non-negativity of the Effective Gravitational Cou- 
pling: 

"0 _ fx^ ^ofxx_ Q 



2«2 



2k^ 



Recall that, for an /(Ar)-theory, the last three require- 
ments are not independent. 



B. Equivalence Of f{X, Y, Z) Gravity With 
Scalar- Tensor Theory 

It has been shown in Ref. [2^ , that the gravity theory 
described by 5 oc / d'^x^/\g\f{X,Y^ Z) is equivalent to 
a multi-scalar-tensor gravity with the scalar fields cou- 
pled to curvature invariants {R, C^, and G)- The equiv- 
alent action, however, is still very complicated and not 
useful to study the particle spectrum of the theory f20|. 
An alternative is to start with the linearized action p2p . 
Actually, by introducing an auxiliary field (/?, the latter 
action can be rewritten in the following form [20| : 



S cx / d^x^\{R- -{Vipf - V{ip) 



2ml 



C^} , (17) 



where 



* Stated in this form Ostrogradski stability implies also, absence 
of spin-2 Weyl ghost propagating modes. 
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and the following scalar field self-interaction potential has 
been introduced: 



3m| 



(18) 



The metric (7^1, has been rescaled as = e'^g^v 

{{S/ifY = g^^V fiipV^ifi). In what follows we consider 
that the ghost propagating mode associated with the 
Weyl squared tensor - last term in the above action - 
decouples from the gravitational spectrum. Otherwise, 
we assume m^^ =0 =4> fy — ~^fz (see the defini- 
tion of in equation (jisp ). The latter requirement is 
automatically satisfied if the invariants Y and Z enter 
in the function f{X,Y,Z) in the following combination: 
Z — AY, as, for instance, in the Gauss-Bonnet invariant G 
[2^ (see also [l^). As long as we wih be concerned here 
with the mass of the scalar mode, this simplification will 
make easier further mathematical handling and physical 
interpretation of the results. 

Looking at the potential ([T51) - where the cosmological 
constant has a non- vanishing contribution - and, since we 
started here with a linearization of the original theory 
and not with the theory itself, it might seem that the 
above equivalence among / d'^x^/\g\f{X,Y, Z) and (fT7)) 
is not as strict as the similar equivalence between a given 
f{R) theory and a dual scalar-tensor theory. At least one 
expects that the mass of the scalar field can be modified 
by the would be ghost excitation, even if it is decoupled 
from the gravitational spectrum p7j . However, as we 
will immediately show, the equivalence is strict. Our 
argument will rest on the computation of the effective 
mass of the scalar field. For this purpose let us to rewrite 
((22)1 as an f{R) = /(Ar)-theory (recall that we have set 
m^^ = 0): ' 



\g\f{X), f{X) = X + -^X'-2A. (19) 



6toq 



It has been shown that, in the weak field limit, looking 
for small spherically symmetric perturbations around de 
Sitter space with constant curvature Xq = Rq, the effec- 
tive mass of the scalar degree of freedom (p = fx in the 
equivalent scalar-tensor theory, is computed though the 
following equation (see, for instance, Ref. (llj): 



Jx ^ ^oJxx 



XX 



of the theory, and on a series expansion around maxi- 
mally symmetric spaces of constant curvature, ii) it is 
based on a particular choice m^^ = 0. The former draw- 
back is common to standard demonstrations of equiva- 
lence between a given /(i?)-theory and its scalar-tensor 
dual, while the latter one is not really a fiag, but a neces- 
sary condition for the absence of Weyl ghosts in general 
/(A:,y,Z)-theories. 



Hence, since for f{X) in ([T^ : 



Critical Comments 



Here we will take a step aside, to discuss the impor- 
tance of a careful and consistent investigation of the pic- 
ture generated by the expansion around maximally sym- 
metric spaces of constant curvature discussed above in 
this section, and the consequent calculation of the masses 
of the excitations associated with the linearized theory. 

At this point the first topic we want to comment on 
is the one related with the lack of coincidence in the ex- 
pressions one finds in the bibliography for the masses 
(squared) of the additional propagating spin-2 and spin- 
degrees of freedom m| and jtiq (the same applies to the 
cosmological constant A emerging after linearizing ([5])). 
In Ref. [2y], for instance, in the definition of the inverse 
mass squared coefficient jtiq ^, there appear terms oc R^ 
so that, when substituted in (IT6t these give contributions 
of the fourth order in R, while the remaining coefficients 
rn^^ and A produce quadratic contributions only. This 
points to an inconsistency of the order of the expansion 
around background spaces of constant curvature i?o con- 
sidered therein. Notice that if the terms oc R^ in the 
expression for TOq ^ in Ref. [26j are removed, the results 
are consistent with the ones reported above. 

Another example is supplied by Ref. 19]. In that 
reference the authors show expressions for the coefficients 
A and m^^ that are consistent with an expansion keeping 
all orders up to i?*, while is consistent with keeping 
only up to R^ . The crucial point in this case is that, if 
in the expansion (|9|10|) one keeps all orders in i?, i. e., 
up to terms ~ i?^ - which is legitimate if one considers 
i? ~ 1 - then, not only the masses mg and will receive 
additional contributions but, at the same time, there will 
necessarily appear in (|16|) additional terms proportional 
to R^ and i?* of the following form: 



FO 
IX 



Xq 
3m2' 



FO 

I XX 



3m§ 



then the effective mass of the scalar degree of freedom 

coincides with the mass of the spin-0 excitation mg com- 
n2 _ ^2 



puted through p^ : ?7T-e// = "^o- 

The above demonstration has several flags. Amongst 
them are the following: i) it rests on the weak field limit 



°^34„2n2 , 044^2 "^54 ^„2 

--^c R + —g ~—gR 



where the coefficients are given by: 



—eg 

(20) 



6 



ai3 






t 023 = fyx 


+ V°ZX, 






^33 


~ fvx ' 


fzx^ 


Oi4 = fyy + 


- »fzY + 






024 


= fvY " 


^Gfzy 


+ ^fzzi O34 


= fyy^ 


- 5/Iy 4 




044 


= fry " 


^^fzy 


+ 4/^z. 054 


-fyy^ 






064 


= fry ~ 


^Vzy 


+ fzz- 






(21) 



about maximally symmetric spaces of constant curvature 
X = = i?o, up to ~ (X* — ^o)^, and keeping up to 
~ Xq. As discussed in section II, the latter is an 

independent requirement having nothing to do with the 
order of the expansion around Xq,Yo, Zq. One obtains 
the following expression to determine the masses of the 
propagating spin-0 and spin-2 excitations (compare with 

Eq. ^y. 



These additional contributions towards (jl6l) are prob- 
lematic and, as long as we know, their inclusion in the 
linearized action has not been investigated in detail yet, 
so that we do not really know how to deal with them. 
Notice, in between, that some expressions in 26] and in 
coincide, while others do not. 

To conclude these critical comments we want to un- 
derline that, in general, it is misleading keeping terms 
with order higher than in Taylor expansion, since, in 
particular, it arises an illusion that one can extrapolate 
the obtained formulas to get to qualitative conclusions 
about the behavior of the masses of the excitations at 
large R-s. As we have shown, this is a wrong procedure 
since the definition of the masses of the excitations is 
given following the study of the linearized action [l^ [l3| 



S(x / d^xVlsl -2A + i? 



6ml 



2ml 



(22) 

while consideration of terms oc R^ ,R^ yields that the 
above linearized action is complemented with the addi- 
tion of the terms in pO)). The actual situation is in fact 
a bit more dramatic; at large R-s the linearized action 
above is replaced by one of the following form: 



where a, /3, 7, 5, and /i are constant parameters. The lat- 
ter action shares nothing in common with the linearized 
action above. 

III. DIRAC-BORN-INFELD LAGRANGIANS 

In this section we shall explore the consequences of ap- 
plying the DBI strategy to /(X, F, Z) theories, from the 
particular viewpoint of their stability properties. In this 
regard it is of prime importance to study the perturba- 
tions of the theory about background spaces of constant 
curvature. In fact it suffices to consider expanding the 
action 



(A + 2/o)(/° -Xo/^^)+Xo(/^)^ 
(A + 2/o)(3/o^ + 2/o+2/|)-3(/0)2 
(A + 2/o)(/° -Xo/°^)+Xo(/i)^ 
(A + 2/o)(/^-f4/0) 



(24) 



In the formal limit A — >■ 00 - unbounded curvature - one 
recovers the known expressions for mg and m| in Eq. 

Notice that, independent on the scale A, the necessary 
condition to avoid the Weyl ghost is the same as before: 
Jy = —4:fz^ so that the DBI modification plays no role in 
avoiding the occurrence of massive, spin-2 ghosts modes 
arising in non-linear theories of the kind f{X, Y, Z). 

The expressions for the masses of the propagating 
linear modes are important since the stability of these 
modes depends crucially on their sign. Hence, in the first 
place, the DBI modification affects the stability prop- 
erties of the gravitational spectrum of the propagating 
excitations. Under the DBI modification of the original 
f{X, Y, Z)-theory, the trace equation ([T4| is replaced by 
the following: 



2A 1 



2fo 
A 




(25) 

Once again, for unbounded curvature A — )■ cxd, equation 
p4| is recovered. The above trace equation is usually an 
algebraic equation to determine the curvature Xq of the 
maximally symmetric background space about which we 
are perturbing the original equations of the theory. 

In the next subsections A and B, we investigate a cou- 
ple of examples taken from the literature to illustrate the 
discussion above. In fact we shall examine only the ab- 
sence of tachyon instability and of the Weyl ghost, which 
are just an aspect of the stability of these theories. Then, 
in subsection C, we discuss in more detail other relevant 
aspects of the stability issue, and apply the results to 
these examples. 



A. Example I 

Take as an example the f{X,Y,Z) theory studied in 
(proposed in Q in a cosmological setting): 



Sdbi oc 




(23) 



fiX,Y,Z)^X 



An+2 



[aX^ + bY + cZy 



(26) 
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where is a parameter with dimensions of mass, and n is 
positive. In the present example the trace condition 
for a maximahy symmetric space of constant curvature 
X = Xq, amounts to 

" ^ ' y'^') 

where we have introduced the parameter k = a + 5/4 + 
c/6. Hence, considering large or small curvatures is not at 
will. Actually, given a fixed set of overall parameters a, 6, 
c, n, and fi, the value of the curvature Rq is fixed through 
([27)) . which is the unique root of the trace equation (fT4)) 
for f{X, Y, Z) given by (PSI) . In this case, for the masses 
of the propagating excitations one obtains, according to 



where A = = Mpi is the maximum curvature scale, 
while a, and (3 are overall constants. Notice that the 
"unmodified" action S (x J d'''x^/\g\f{X,Y, Z) coincides 
with the Einstcin-Hilbcrt one, since the Gauss-Bonnet 
term X^ — 4Y + Z contributes a total divergence that can 
be safely erased from the action. Hence, in the linearized 
approximation, there is only one massless spin-2 propa- 
gating excitation (properly the graviton), while the grav- 
itational spectrum of the theory depicted by the action 
([29l) . in the linearized limit, consists of a massless spin-2 
propagating degree of freedom, and of a spin-0 (massive) 
excitation. Otherwise, it is properly a scalar-tensor the- 
ory of gravity. In this case the DBI procedure evidently 
affects the "unmodified" f{X,Y,Z) theory, even at the 
linearized level. The mass of the scalar degree of freedom 
in the linear approximation is given by (see Eq. (j24p ): 



777-2 = p 



n+1 

nl 
n+l 

nk 



(fc2 -I- 2na^)^2 ^ 
fc2 -I- 2na'^ " 



b + ic 



(28) 



where 



and 



I = {3a + b + c)k - 6(n -|- 1)0^, 



2(n+ 1) 



l/(2n+l) 



It is curiously enough that, in the present case, for the 
corresponding DBI modified theory (P^ . the expansion 
around spaces of constant curvature Xq yields the same 
results for the masses of the spin-0 and spin-2 excitations 
as in Eq. (gS]).^ This means that, for the theory (1^ . 
up the given approximation, the stability of the pertur- 
bations is not affected by the DBI modification 



B. Example II 

To illustrate how the DBI modification does evidently 
affect the dynamics, we choose a theory that fulfills the 
requirements to be phenomenologically viable, listed at 
the beginning of the section. It is based on the following 
DBI-type action, which represents a minimal deformation 
(a trivial change of signs only) of the one in j26|] : 



5 cx 1 d^x^l^ [^l + 2K^f{X,Y,Z)-l) , 

f{X, Y, Z) = aX - K^(3{X^ -AY + Z). (29) 



^ To be consistent with the order of the approximation, one has to 
^•^^"^ and higher. 



I + SaK^Xo 

3K2(4/3K2Xo-a)^ 



(30) 



where wc have dropped terms of order cx Xg, k^Xq and 
higher. The value of the background curvature is a real 
root of the trace equation ([25]) . Notice that in the un- 
bounded curvature limit — > (A — > oo), the mass of 
the scalar perturbation blows up and is negative, so that 
this limit, if it exist in the present theory, is largely un- 
stable. Besides, in order for the theory (1^^ to be free of 
tachyon instability (negative m^-s), unless a is negative, 
the value of the background curvature is bounded from 
below 

Xo > 



On the other hand, the constant a can not be negative 
because, in the low curvature limit, it modifies the sign 
of the the effective Newton's constant G — {8n)^^ Mp^ . 
This demonstrates that the theory given by (|29)) does 
admit stable perturbations around maximally symmet- 
ric spaces of constant curvature (otherwise, it does ad- 
mit stable de Sitter solutions), only for unnaturally small 
values of the constant a: a ^ Pk^Hq = PMpfH^ {H^ 



neglect terms oc ~ 



is the current value of the Hubble parameter), so that 
we are faced with a very serious problem of fine tunning. 
Although possible, this is a very unlikely scenario. Of 
course, the situation here is that the stability of the the- 
ory is very sensitive to the signs in ()1|23|) . If one chooses 
an alternative deformation a la Dirac-Born-Infeld as it 
follows (compare with (jlj or, equivalently, with (|23p and 
note the change of signs) : 



/(X, y, z) ^ A (i - v/1 - 2/(x, y, z)/\) (31) 

then, after expanding the corresponding action around a 
maximally symmetric space of constant curvature X ~ 
Xq, one is led to the following expression to determine 
the mass of the spin-0 excitation (compare with (IMl) ): 
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(A-2/o)(3/o^ + 2/o+2/i)+3(/0)2 



(32) 



In the theory with f{X, Y, Z) given by (|29l). one obtains 
(compare with Eq. ([5(11)): 



1 - ian'^XQ 
3K2(a-4/3K2Xo) 



(33) 



where, as before, to be consistent with the approxi- 
mation undertaken, we have dropped terms of order 
oc XqjK^XqjK^Xq and higher. The tachyon instabihty 
is absent if either: i) Xq < l/^K'^a and, at the same 
time, Xo < a/4:K^j3, or ii) Xq > l/^K^a and, at the 
same time, Xo > a/AK^p. The latter possibihty is very 
unlikely to occur since, to obtain a stable scalar per- 
turbation, the constant a has to be unnaturally small 
a ~ jSMpfH^ (see similar discussion above). The for- 
mer possibility is indeed appropriate since it is compat- 
ible with a maximum curvature, as it is expected for 
DBI models. For > /3, the limiting (larger avail- 
able) curvature is X™°^ = l/Sn^a, while for /3 > 
it is X^°-'^ = a/4:K?l3. So, perhaps, the alternative DBI 
strategy ([31]) is more attractive than the first one (P) in 
the present case. 



C. DBI Modification of Stability 

Here, we will summarize the way the DBI strategy 
modifies the stability requirements listed in subsection 
A of section II for arbitrary functions f{X,Y,Z). For 
this purpose we write the two possible choices of the DBI 
modification in the following general form: 



/^eA(Vl + 2e//A-l), 



(34) 



where e = ±1 (e^ — 1). Written in this form, equation 
([M)) comprises both: procedure given by Eq. ([T]) if one 
chooses e = and the alternative procedure given by 
Eq. (|31l) if one chooses e — ~1. The following relation- 
ships are useful to study the considered modifications of 
the stability: 



fx^rfx, /y^r/y, fz^rfz, 



fxx r i^fxx - ^ fx) ) 
where the DBI "boost" is defined as 

r = r(x,y,z)= (i + y./ 



-1/2 



Hence, after applying the DBI deformation to quadratic 
modifications of Einstein-Hilbert theory, the following 
modifications of the stability requirements occur (com- 
pare with similar requirements in subsection A of section 
II): 



Ostrogradski Stability: The same as before, i. e., 
it is unaffected by the DBI procedure 



/o +4/0 =0 



Ricci Stability: 



Jx 



XX 



j^ifxf > 



• Absence of Tachyon Instability: 

2 fx - ^afxx + J^o^lifx)^ 



mx + 2./^ + 2/0 



fr2(/o)2 



> 



Non-negativity of the Gravitational Coupling: 



Jx 



Xof"xx + lxorl{fxf>o 



Additionally it is required that the DBI boost be a real 
quantity, i. e.. 



Tea? 



2ef + A > 



Notice that, in the formal limit A — >■ oo, the above re- 
quirements on the stability coincide with the ones in sec- 
tion II A. 

Let us to check the theory ([26]) [l^, regarding other, 
perhaps more restrictive, stability criteria. For instance, 
since in order to avoid propagating Weyl ghosts b = 
—4c => k = a — 5c/6, then 

_ 2na[(2n+l)a + 5c/6]//"+2 
fxx — 



<o. 



therefore, for the space of parameters allowed by the ob- 
servations, the theory of Ref. [l^ is Ricci unstable, so 
that it develops ghosts. This argument might be enough 
to rule out this theory. Next see what happens if we 
modify it a la Dirac-Born-Infeld. Would it be yet Ricci 
unstable? After applying the DBI deformation proce- 
dure (p4|) . the above condition on fxx translates into 
the following bound: 



X^" [Xk + ea{n + l)Xo] < 



z[6{n+ l)a + 5c] 
3fc" 



4n+2 



where k = (2n + l)a + 5c/6. For small Xq ^ A, we have 
e[6(7i + l)a- 



Xri 



< 



4n+2 



3kk-^X 

which is consistent only for DBI embedding with e — +1. 
Notice that in the formal limit A — oo the above bound 
can not be satisfied and we get again a Ricci unstable 
theory. For the opposite embedding, i. e., e = — 1, the 
only way to achieve Ricci stability, is that Xq > Xk/{n + 
l)a, leading to the following constraint: 



X, 



-[(n + l)aXo-Afc]>^^il±^/"- 



In this case, A is not an upper curvature bound any more, 
and the original motivation of the DBI strategy is lost. 

We see that the DBI modification strategy indeed mod- 
ifies the Ricci stability criterion as well, making a origi- 
nally Ricci unstable theory, a stable one. 
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IV. DBI DEFORMATION STRATEGY APPLIED 
TO f{R) THEORIES 

We have seen how the apphcation of the DBI proce- 
dure - both variants: ((J), and (PT|) - actuahy modifies 
the stability properties of f{X,Y,Z) gravity theories^. 
Even the gravitational spectrum of the theory is altered 
if one follows this strategy with a class of theories given 
by f{X, Y, Z) = X+aQ (these include the particular case 
a — 0). Here we will be concerned with possible modi- 
fications of the stability of de Sitter backgrounds in the 
simpler /(X)(or f{R)) gravity. In what follows, just for 
homogeneity of writing, we keep the variable X instead 
of i?, although both will be used interchangeably. Be- 
sides, for generality of the discussion we will consider the 
DBI deformation procedure given in Eq. The field 

equations that can be obtained from an action principle, 
the action being: 



5 = 



2^2 



(35) 



-2A = ;/o_^o-^o + X'^fo 1'^ 



^ofo 



(40) 



The gravitational spectrum consists of a massless spin- 
2 graviton, and a massive spin-0 propagating degree of 
freedom, with mass squared uiq. In this case, the trace 
equation (1551) translates into an algebraic equation to de- 
termine Xq: 



Xofo — 2/o . 



(41) 



As it has been already explained, requirements of 
Ricci stability, absence of tachyon instability, and non- 
negativity of the effective gravitational coupling are not 
independent requirements in the present case, so that one 
has to check only two of them, for instance, Ricci stabil- 
ity, and non-negativity of the gravitational coupling: 



where S'(m) is the action of matter, and -0 - the collective 
name for the matter degrees of freedom, can be written 
in the following form: 



nx) g^^^k' 



(36) 



where the comma denotes derivative with respect to the 
curvature X, and we have introduced the following defi- 
nition of the effective curvature tensor: 



K^T, 



^ f{X)-Xf'{X) 



-V^V,/'(X) - g^,Uf'{X) . (37) 



The trace of Eq. p6p amounts to an additional constraint 
on X: 



Xf{X)-2f{X) + 3nf'{X) = 



(38) 



If one expands the action ((35|) - considering vacuum 
background, i. e., S'(m) = - around maximally sym- 
metric spaces of constant curvature X = Xq, just as it 
has been done in section III, then one obtains: 



/">0, /^-Xo/^'>0 



(42) 



The same linearization around vacuum backgrounds 
with constant curvature, when applied to f{X) gravity 
which has suffered further DBI deformation given by Eq. 
([34|) . yields to the following expression for the mass of the 
spin-0 propagating excitation (compare with Eq. (j40l) ): 



2 _ /o ^ ^ofo + ■|rQXo(/o)^ 



where 



3/^'-fr2(/^)2 



v/l + 2e/o/A 



(43) 



while the trace equation (|4T|) translates into the following 
equation: 



TlXofl, - 2eA(l - To) . 



(44) 



Therefore, the stability criteria (|42|) translate into the 
following requirements: 



Socao [ d^xJU f -2A + X + -^X 
J V 6"io 

where (/o = f{Xo), etc.) ao = /o - ^o/o , and 



(39) 



® This is the usual case in DBI modifications, see[35| for an exam- 
ple of a DBI-modified cosmological scalar field. 



/o 



Xof[; + jTlXoif[,r>0 



(45) 



In what follows we will explore the consequences for 
stability of - constant curvature - vacuum background 
linearization, of applying the DBI procedure given by 
Eq. ([M)) to f{X) theories of gravity. We plan do that by 
working out two examples taken from the bibliography. 
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A. Example I 



B. Example II 



We start the study of the possible modification of the 
stabihty of (constant curvature) vacuum background hn- 
earization of f{X) theories - which have been further 
modified according to the DBI deformation strategy -, 
by choosing the model of reference [6^ , that is given by: 



Our next example also reflects the consequences of ap- 
plying the DBI deformation procedure (|34l) . Let us con- 
sider the following f{X) gravity theory [32i] : 



/(X) = X-(l-n)^2(i-n)^« 



(47) 



f{X) = X 



2(n+l) 



X^^ 



(46) 



where /i is a suitably chosen parameter. A check of the 
stability criteria (|^^ shows that 



f" 



n(n + 1)^2(7^ 



X, 



n+2 



is negative, which means that the model suffers from 
Ricci scalar instability [28l - l30j (see also [ill, IHl)- 
Nonetheless, 



/o ^ ^ofd - 1 



n(n + 2)^2(n+i) 



XJ; 



> 



so that positivity of the effective gravitational coupling is 
guaranteed. As a consequence (see the definition of the 
spin-0 mode mass squared in Eq. (I40p ). the model suffers 
also from the scalar tachyon instability. As explained 
before, the Ricci instability is enough to rule out the 
model since this type of instability develops very quickly. 

What one should expect from applying the DBI strat- 
egy in this case? The answer can be based on the anal- 
ysis of the DBI modification of the stability criteria in 
Eq. According to these, the Ricci stability bound 

is given, in the present case, by the following constraint: 



-^0 



n(n-t- l)/i2("+i) 

br small cu 
straint simplifies 



X^ 



2{n+l) 



xr' 



2 > 0. 



For small curvature X^^'^ <C A/i^'"^^-', the above con- 



eA 

'x^ 



1 



2{n+l) 



X, 



n+1 



2 > 



so that for the DBI embedding with e = 



Xq~^^ < ^^("+1^2, this bound is always satisfied 



1, whenever 
Be- 
sides, if in order to simplify the analysis, one assumes that 
Xq ^ /i^, then, the criterion requiring non- negativity 
of the effective gravitational coupling upon linearization, 
can be written in the form of the following bound: 



A 



-4\ ^2(n+l) 
^0 



> 



which for the e = — 1 embedding is always satisfied. 
Therefore, the DBI deformed theory is not only Ricci 
stable, but also, it is free of scalar tachyon instability. 



where, to be compatible with observations < n < 0.25, 
and /Li is a sufficiently small parameter (^ ^ lO^^^eV) 
[33| . If one linearizes the corresponding action 5* cx 
/ df^x^/\g\f{X), then one finds that the mass of the spin- 
excitation Eq. (HU)) . in the present case is given by: 



Too 



{Xofi- 



-2\l-n 



n{n~ l)(n- 2) 



f , (43, 



so that, for positive n-s within the range where the 
model is compatible with the observations (see above) , a 
tachyon instability develops (toq < 0). Notwithstanding, 
since in this case, for the mentioned range of parameters 
compatible with observations 



= n{n - lYfi 



2 2(1-1- 



then the model is Ricci stable. One might think that 
a large enough time to develop the tachyon instability 
could make the theory compatible with the cosmological 
dynamics. However, the very tricky situation with this 
model, comes, precisely, from the fact that Ricci stabil- 
ity, absence of tachyon instability and non-negativity of 
the effectve gravitational coupling are not independent 
requirements. In this particular case (Ricci stable, but 
with tachyon instability), the effective gravitational cou- 
pling is negative, which is catastrophic for this model, 
and it has to be ruled out. 

Would replacement of f{X) by eX{.yi + 2ef{X)/ X- 1) 
save the model? Let us start by checking Ricci stability. 
In this case, since /i^ is small enough, one might consider 
two important limiting situations to simplify the analysis. 
Assume first that Xq ^ /i^, so that the Ricci stability 
bound in (l45l) can be written as 



which is satisfied only for the DBI embedding with e = 
-|-1. Now assume a different limit, Xq ~ In this 
case, the above bound may be written in the form of the 
following constraint: 



en(l - nfX - pXo > ^ Xq < 



en{l 



-X 



P 



where the constant p = I + n{l — n)^ — 2n{l — n)(2 — n) 
is positive for the range of parameters allowed by the 
observations. This means, once again, that the only DBI 
embedding that is enable to save the Ricci stability is 
the one with e = -1-1. The next step is to check for 
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non-negativity of the effective gravitational coupling. In 
this case, for the limiting situation when Xq ^ /i^, the 
criterion of non-negativity of the effective gravitational 
coupling in (|45p amounts to the following inequality: 

vr^ < ^ 

° - 2(l-n)(l + 3n-?i2) ' 

which is obeyed whenever e = +1 as it should be ac- 
cording to the previous results on Ricci stability. If one 
considers, instead, the limiting situation when Xq ~ /i^, 
one obtains that the effective gravitational coupling is 
non-negative if 

Xo < , 



where the constants fc = 1 — n{l — n)(2 — n), and I = 
2(1 — n){l + 3rt — rt^) — n(l — 71)^(4 — n), are positive 
for the range of parameters allowed by the observations. 
Hence, once again, the bound is satisfied if e = -1-1. 

We have seen that, in general, the DBI strategy to 
modify f{X) theories can help avoiding either ghosts or 
the tachyon instability or, it can even help avoiding both 
simultaneously. For Lagrangians where the quadratic 
contributions come with inverse powers of the curvature, 
the correct DBI embedding is the one with e = — 1, while 
for quadratic contributions proportional to positive pow- 
ers of the curvature, the correct embedding is the one 
with e = +1. In the next section we will study the modi- 
fications to the cosmic dynamics caused by applying the 
DBI strategy to f{X, ...)//(X)-theories. 



COSMOLOGICAL CONSEQUENCES OF 
APPLYING THE DBI STRATEGY 



In this section we will explore the impact of applying 
the DBI procedure - both possibilities: (|ll23p . and ((3T|) 
- on the cosmological dynamics of f{X, ...)/ f{X) gravity 
theories. Due to their simplicity, we focus the discussion 
in /(i?) gravity theories exclusively. For the study of 
the impact of the DBI strategy applied to f{X, ...) the- 
ories, we refer the reader to the paper [34;!, where this 
issue has been investigated in detail, for the case when 
the above deformation procedure is applied to theories of 
the kind f{X,0). We choose the Friedmann- Robertson- 
Walker (FRW) metric as a local description of spacetime 
at cosmological scales (for simplicity, without loss of gen- 
erality, we consider FRW spaces with spatial sections of 
constant curvature):^ 



(49) 



Inserting this choice of the metric into the field equation 
(|36l) . and assuming a perfect fluid stress-energy tensor 

where is the fourth-velocity of an observer co-moving 
with the fluid, p„j and pm being the energy density and 
the pressure of the fluid, respectively, to describe the 
macroscopic behavior of matter, one obtains the follow- 
ing cosmological equations: 



3^2 



/ 



I {Prn ~t~ Pcur) 



2H + = - — {prn + Pcur) , (50) 

where the effective (parametric) energy density pcur and 
pressure Pcur of the effective "curvature" fluid, are given 
by: 



1 



Hcur — 2 
hi 

1 

Pcur 



-{Xf~f)~3Hf 



2Hf' + f'--{Xf'-f) 



(51) 



respectively. As it is clear from going to the limit Pm — ^ 
in the first equation in (|50p - properly the Friedmann 
equation -, the effective energy density pcur can not be 
negative [11|. It can be defined also an effective equa- 
tion of state (EOS) parameter Wcur = Pcur/ Pcur in the 
following way: 



Wc 



2(/' - H,h 
Xf'-f- 6Hf' 



(52) 



from which one may see that the vacuum value Wcur = 
— 1 is attained whenever 

/' = Hf' . 

de Sitter solutions trivially satisfy this equality, but there 
is a class of non-dc Sitter solutions for which /' cx a{t) 
which also satisfy the above equation. 

Under the replacement (l34l) . /', /', and /', transform 
like: 



For a discussion of the subtleties associated with the validity /' — > F ■( /' — f'^ ( X + i—X — 3^^— /'X^ 

of the Ehlers-Geren-Sachs theorem |3h| on the identification of 



FRW spacetime with our universe, as a consequence of the high 
degree of isotropy of the cosmic microwave radiation, we refer 
the reader to and linked references therein. 




F/' , (53) 
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respectively, where F = (1 + 2e//A)~^/^. Hence, the 
effective energy density and pressure of the "curvature 
fluid" (Eq. (fSTjl ). wiU be replaced by very complicated 
expressions. For that reason, to simplify the analysis, we 
will operate in a different, yet equivalent, way: we keep 
equations ([50| . ([5T|) . ([52]), valid for some fe{X), such that 



MX) = eA(Vl + 2e/(X)/A - 1) , 

where f{X) is the original, unchanged function, while 
fe{X) will be the DBI modified one. Otherwise, to look 
for the impact of the DBI modification strategy on the 
cosmic dynamics of f{X) models, one has just to replace 
f{X) — ^ fe(X) in the mentioned equations. In the next 
subsection we will investigate the dynamics of given f{X) 
models, to uncover the way the DBI procedure operates 
to modify it. In what follows, for simplicity, we will con- 
sider pressureless dust matter as the background fiuid, i. 



e., Pr, 



0. 



A. Dynamical Systems Study 

Now, since finding exact solutions of the equations ((50)) 
is in general a very difficult task, we will rely on the 
dynamical systems tools to investigate the asymptotic 
structure of the f{X) models of interest, instead. To 
this end we will apply the concise recipes given in the 
appendix (section VIII). The goal will be to write the 
system of cosmological equations in the form of an au- 
tonomous system of ODE (as described in the appendix) , 
so that one could associate such important dynamical 
systems concepts as past and future attractors (also sad- 
dle equilibrium points), with dynamical configurations - 
solutions - of the models. This is a powerful approach to 
uncover the most generic classes of solutions that are al- 
lowed by them. In order to build an autonomous system 
of ordinarydifferential equations (ODE) out of (|50l5ip . 
following [S^l, we introduce the following dimensionless 
variables: 



X 

6IP 



H 



y 



Hf 



(54) 



In terms of these variables the Friedmann equation (first 
equation in (|50p ) can be written in the form of the fol- 
lowing constraint: 



^7D, 



2 

1^ Pm 



= 1 



6iJ2/' 



(55) 



where we have conveniently defined an "effective" di- 
mensionless matter energy density parameter O^, which, 
thanks to the /' function entering in its definition, might 
be, in principle, any sign, i. e., either positive or nega- 
tive, without entering in confiict with standard physical 
requirements. Particular properties of a given model are 



encoded in the last term in the RHS of Eq. (|55|) . In 
fact, not all of the variables of the phase space needed 
to describe the dynamics of a given model have been de- 
fined yet. This requires knowledge of the concrete model. 
Other general, useful expressions are 



/' 



2ix - 2) + y + eClr, 



= -1 - 



2(a;-2)-H60^ 



///' 



y 



(56) 



The phase space variables x, y obey the following au- 
tonomous system of ODE: 



dX/dr 
6iJ2 



- 2x{x - 2), 



dx 
d^ 

^^{x~2){2-y)+y 



y 



(57) 



where, as customary, we have introduced the time- 
ordering variable r = In o - basically the number of e- 
foldings. Depending on the concrete model, new vari- 
ables have to be introduced and, as a consequence, the 
above equations have to be complemented with the addi- 
tion of new ones (one per each new phase space variable) . 
The particular properties of a given model enter in equa- 
tions ([57)) through the terms dX/GH^dr and (l^ (see Eq. 
([55)) V Now we are in position to study a concrete model. 



B. Model f{X) = fiX" 

Although, in general (n 7^ 1), it does not contain 
Einstein-Hilbert gravity as the low curvature limit (as 
it should be), the above model is simple enough, and has 
been formerly studied, for instance in jS^j (IJMPD). The 
related model f{X) = X + jiX"' - also studied in [sl] 
(PRL) - is a bit more realistic, yet much more compli- 
cated from the point of view of its asymptotic proper- 
ties. However, the real complexity arises when one tries 
to study the dynamics of its DBI deformation. This will 
deserve a separate publication. That is the reason why 
we have not included this toy model as an example in 
this subsection. 

In the model f{X) — /iX", /i and n are the overall free 
parameters. In this case, since ///' = X/n^ hence 



y ■ 



Besides, since 



then 



/' = n(n- l)/iX"-2x 



dX ( 1 



dr V 1 — n 
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Using the above relationships, one can write ([57]) in an 
exphcit form: 



dx 
dT 



1 - n 



xy — 2x{x — 2), 



= 61 



+ {x-2){2-y) 



-5y + y\ (58) 

The above is a closed system of ODE, and one does not 
need more variables to describe the phase space dynam- 
ics, otherwise, the autonomous system (|58|) is defined 
within the following (open) 2-dimensional phase space:^ 

^^{{x,y):x>Q, n,n{x,y)>0}. 

The latter bound is justified only if one assumes the pa- 
rameters fi and n to be both positive magnitudes, so 
that /' is always non-negative. Other two magnitudes 
of immediate physical (and observational) meaning are, 
the deceleration parameter q — 1 — x, and the effective 
equation of state parameter of the "curvature fluid" 

2{x-2) + 6n„, 



3(^)x + 3y 



The critical points Pi — (xi, yi) of the autonomous sys- 
tem of ODE ([55|. together with their relevant properties, 
are listed below: 

1. Curvature-dominated equilibrium point Prad — 
(0, 1), associated with decelerated peace of expan- 
sion. The "curvature fiuid" mimics radiation. This 
critical point is characterized by the following val- 
ues of the relevant parameters: 



0, Wc 



1 



q=l. 



The eigenvalues of the corresponding linearization 
matrix are (see the appendix) 



Ai 



-1, Ao 



4n - 



n-l 



For 1 < 71 < 5/4 the point Prad is a future at- 
tractor (not adequate to the present cosniological 
paradigm) . Otherwise, Prad is a saddle equilibrium 
point. 



2. The critical point 



p — 



2n- 3 6(n- 1) 



is associated with matter-curvature scaling. It is 
characterized by the following values of the relevant 
physical parameters 



q ■ 



3- lln + 7n2 



3(2n- 1) 
8n-3 



while the eigenvalues of the corresponding lin- 
earization matrix can be written, in compact form, 
as it follows: 



Ai, 



3n2 - 6n -1- 3 ± 

2n{n- 1) ' 



where s ee 81 -t- llQn^ - 420n + 121n'' - 49271^. The 
above point exists, i. e., it belongs in the phase 
space, whenever 0.35 < n < 1.3. In other words, 
for n-s within the latter narrow interval, Clm > 0. 
Besides, for the allowed range of the parameter n, 
q > always, so that Pm~c is associated with de- 
celerated expansion. For 0.35 < n < 0.38, and 
1 < n < 1.22, the equilibrium point Pm-c is a sad- 
dle critical point, while for 0.38 < ri < 1, it is a 
stable spiral (future attractor). In the latter case, 
since <? > 0, the model is not suitable to accommo- 
date the present cosmological paradigm. 

3. The curvature-dominated critical point (f2„i — 0): 



Pc-d = 



n(4n-5) 2(n-2) 
{2n - l){n - 1)' 2n - 1 

is characterized by 



'^r.ii.r 



6n — 7n — 1 
3{2n^ - 3n-f 1)' 



2n^ -2n-l 
2n2 - 3n -I- 1 ' 



and the following eigenvalues of the Jacobian ma- 
trix are obtained: 



Ai 



2(7n2 - lln + 3) 
2n2 - 3n + 1 



A, = 



An -5 
n-l' 



Here we will consider only positive scalar curvature. 



The point Pc-d can be associated with accelerated 
expansion for the following ranges of the free pa- 
rameter n: 

1 , l + \/3 

- < n < 1, < 7T, < oo , 

2 2 

i. e., practically for all of the real segment (but for a 
very narrow interval). Meanwhile, the above eigen- 
values are both simultaneously negative for the n- 
range: 5/4 < n < cx). In other words, the (al- 
most every where within the n-interval) inflation- 
ary equilibrium point Pc-d is stable - a future at- 
tractor - but for the narrow interval < n < 5/4. 
For the particular value n = 2, since x — 2 =^ 
H = 0, this equilibrium point is associated with a 
de Sitter solution. 
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Notice that, given a value of the variable y = yi in a 
critical point Pi = {xi,yi), then, due to the definition of 
this variable, one can write the curvature scalar X as a 
function of the scale factor 

X{a) = Ma-^'/("-i) +Xo , 

where M and Xq are integration constants. Hence, for 
instance, for the critical point Pc~d, 

2(n-2) 

X{a) cx o(2"-i)("-i) . 
This means that, for 

n< -, 1 < n<2 , 
2 

the scalar curvature decreases with the expansion of the 
universe, while, for 

i<n<l, n> 2 , 

the curvature unboundedly grows as the expansion pro- 
ceeds, thus mimicking phantom behavior as a mere cur- 
vature effect. 

In general, the above results coincide with those in 
[3^ . but for the point associated with the so called by 
the authors "0 matter dominated epoch" (their point 
C), which, in our study corresponds to the equilibrium 
point Prad, where the non-linear curvature effects mimic 
radiation. Besides, several expressions for the matter- 
scaling solution 



32| 



- do not coincide 
321 chose a differ- 



solution B in 

with ours. Even if the authors of Ref. 

ent phase space variable X2 = — //6/'iJ^ (instead of our 
X = X/6H'^), nonetheless, the results should coincide. It 
seems to us that the variable X2 in [32|] is a bad choice 
after all. Actually, in the first stages of the study pub- 
lished in Ref. [3J| , we started by using X2 as one of the 
phase space variables. Unfortunately, a problem arose: 
the limit in which the unmodified f{X, Y, Z) is recovered 
from its DBI modification / {X, Y, Z) was not obtained! 
We were then forced to renounce to the variable X2 ■ This 
was, precisely, the motivation to use a different variable 
in fSJ], and in the present investigation as well. We rec- 
ognize that this is a tricky situation, and that an inde- 
pendent investigation of the nature of this inconsistency 
is necessary, but this is behind the scope of the present 
paper. 



C. DBI Modification of the Dynamics 

In the present subsection we aim at investigating the 
way the above dynamical systems picture is affected by 
the Dirac-Born-Infeld strategy. Recall that, in order to 
consider the DBI modifications to the model in the above 
subsection {f{X) = fJ,X"-), in the equations (I54II57P one 
has to replace 



f{X) ^ f{X) 




so that, for instance (compare with equations (|54II57I) ) 



nr. 



=l-x-y 



2{x-2) + y + 6n,n, 



U)r,ii,r — 1 



2{x - 2) + 6nr, 



x-^+y 



(59) 



etc. 



The first consequence of the DBI modification of 
f{X) = /xX", is that the system of ODE ^ is not 
a closed one, so that, in addition to the variables 



X 



y ■ 



one has to consider the new phase space variable 

eA 



z = 



Hence, for instance 



z(/3-l)/3 



where (3 — f3{x,y,z) = + 2x"'/z. It is also verified 
that 



dX 



z + 2a;" 



(n - l)z +{n~ 2)x''' 



xy 



After appropriate algebraic manipulations one obtains 
the following autonomous system of ODE for the DBI 
modified version of the model f{X) — fiX": 



dx 

d^ 

dy 

dT 
dz 

d^ " 
where now 



z + 22:" 



{n - l)z +{n- 2)a;" 
{x-2){2-y) + y + y^ +6n 

-2n(x- 2)z, 



xy — 2x{x — 2), 



(60) 



nm = l- x-y + 



/3(/3-l)z 



— 1 



The corresponding (open) 3-dimensional phase space 
where to look for critical points of the above system of 
ODE, can be defined as 

* = {{x, y,z): x>0, Cl,n > 0, 2x" + z > 0} . 

The latter bound comes from requiring that > 0, while 
the former one ilm = Pm/QH^ f > 0, is justified only if 
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TABLE I: Results of the numerical study of the eigenvalues of the linearization matrix corresponding to the equilibrium point 
Pt,„ in subsection C of section V. 



71 




y 


z 


^1 




^3 


2 


2 





-1.44 


-8.64 


2.54 


-2.89 


2.5 


)) 




-1.04 


-3.6 + 4.5i 


-3.6-4.5j 


-1.8 


5 


i5 


■>7 


25.5 


-5.76 


-3.95 


0.7 


10 


)) 


1? 


3144 


-6 


-4 


1 



TABLE II: Results of the numerical investigation corresponding to the equilibrium point P, in subsection C of section V. The 
values of the free parameter n are approximately concentrated within the range 2 < n < 3, since, for other values, the effective 
dimensionless energy density parameter Clm is negative, meaning that the corresponding point does not belong in the phase 
space. The value n = 1 is an isolated one. For n = 3, P, is a non-hyperbolic critical point. 



n 


X 


y z 


Ai 




A3 






q 


1 





1 


6 


-1 


4 





1/3 


1 


2.02 


1.17 


0.02 " 


-3.2 + 13.5i 


-3.2 - 13. 5j 


3.33 


0.27 


-0.98 


-0.2 


2.05 


1.23 


0.04 " 


-3.3 + 9.62i 


-3.3-9.62i 


3.16 


0.24 


-0.97 


-0.23 


2.1 


1.33 


0.07 " 


-3.4 + 6.9j 


-3.4 - 6.9i 


2.8 


0.2 


-0.95 


-0.32 


2.3 


1.61 


0.11 " 


-3.8 + 3.86i 


-3.8 - 3.86i 


1.8 


0.1 


-0.93 


-0.6 


2.5 


1.8 


0.1 


-4.1 + 2. 7j 


-4.1 - 2.7i 


1.1 


0.05 


-0.95 


-0.8 


2.7 


1.9 


0.07 " 


-4.3 + 2.01i 


-4.3 - 2.01i 


0.53 


0.02 


-0.97 


-0.9 


2.9 


1.97 


0.24 " 


-4.4+ 1.55i 


-4.4 - 1.55i 


0.15 


0.005 


-0.99 


-0.97 


3 


2 





(-9 + V7i)/2 


(-9 - ^/7^)/2 








-1 


-1 



the free parameters (/it, n), are both positive magnitudes. 
In fact, since /' = P~^f', then, as long as /' > (which 
is true only for positive {fJ.,n) as noted above), hence /' 
is also non-negative. 

Notice that all of the equations and expressions in the 
unmodified case f{X) = iiX" - studied in the former 
subsection - are recovered from the ones above, in the 
limit z — >■ oo. Actually, in this limit: 

dX xy 
hm — = , 

z^oc dT n — 1 

while, since, up to second order in Taylor expansion 
lim /3 = 1 + — + 0{z^^) , 

2-)-00 Z 

then 

z-foo QH^f n \ z ) n 

The relevant equilibrium points of the autonomous sys- 
tem of ODE (pn|) are listed below, and their main prop- 
erties are also shown. 



1. The infiationary, curvature-dominated, de Sitter 
equilibrium point in the phase space \1/ (see the 
definition above), 

P,5 = (2,-l,-2"+i) , 

is characterized by the following values of the pa- 
rameters of observational interest: 

The eigenvalues of the Jacobian matrix correspond- 
ing to Pis are 

Ai - -7, A2 = -4, A3 = -2 , 

so that it is always the late-time attractor, inde- 
pendent of the value of the free parameter n. In 
this case the effective "curvature fiuid" behaves as 
vacuum energy - properly, as a cosmological con- 
stant. 
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2. A twin - also inflationary - de Sitter state in the 
phase space, corresponds to the point 



Ptra = 2,0, 



2"(n-4) + fc + l 



where k = Vl + n2"+i. As for its twin-solution 
above, for this equilibrium point one obtains that 



nr. 



0, Wc 



-1, q 



-1 



This time the magnitudes of the eigenvalues are 
given by extremely huge (and complex) expres- 
sions, so that, in their place, we decided to present 
the results of a numerical study for several values of 
the parameter n instead. These results are shown 
in Tab. HI As seen, Pt^ can be either a saddle criti- 
cal point, or an stable spiral. In the latter case Pft„ 
rivals with its twin-solution Pds'- two latc-timc at- 
tractors co-exist, a feature also found in the study 
presented in SJj. In the former case, due to its 
transient character, the solution Ptw could serve 
as an alternative explanation to primordial (early- 
time) inflation. 

3. There is another equilibrium point 

P* = (s,y,o) , 

where x = x{n) and y — y{n) are complicated 
functions of the parameter n. The corresponding 
expressions for n„i — flmin), Wcur = Wcurin) and 
q = q{n), are also bizarre huge expressions of the 
argument n. That is the reason why we present 
only a numerical investigation of the properties of 
this equilibrium point (the results are displayed in 
Tab. |ll|. In general, whenever it exists, is a 
spiral saddle point, but for n = 1 (see the dis- 
cussion below), and for n = 3, where it is a non- 
hyperbolic point. For the particular, isolated, value 
of the free parameter n = 1, one has P* = (0, 1, 0), 
andOm = (curvature-dominated phase) . Besides, 
since Wcur = 1/3, the curvature behaves as a radi- 
ation fluid, favoring decelerated expansion {q = 1). 
This time, since the eigenvalues of the Jacobian 
matrix are of different signs: 



Ai — 6, A2 



-1, A3 



then, this phase represents a saddle equilibrium 
point in the phase space. For other values n± 5n 
in the vicinity of n = 1, since f2„ is negative defi- 
nite {Vtn < 0), then the point P, does not belong in 
the phase space 5' (see the definition above). The 
delicate situation with this equilibrium point is as- 
sociated with the fact that, since z = 0, then, either 
/J, — >■ 00, or — >■ 00. While the former case is not 
of physical interest, the latter one leads to states of 
unboundedly large curvature (perhaps a singular- 
ity). As long saddle point, which means. 



in turn, that the corresponding state can be only 
asymptotically approached by the system, this is 
not catastrophic. 

In general, since at a given equilibrium point P = 
(x, y, z), y = —f'/Hf, then, one can write the following 
integral in quadratures: 



dX l{n- 1)A -f e{n - 2)^iX''- 



A -f 2e^X" 



= -y In a + Co , 



where Co is an integration constant. Once the above inte- 
gral is computed one gets a ~ a(X). Hence, by inverting 
the latter function (whenever the inverse exists) , one can 
write the scalar curvature as a function of the scale factor 
X ^ X(a). 

Worth noticing that, while for the unmodified /(X) = 
/iX"-theory the case with n = 1 is just Einstein- 
Hilbert gravity with only one critical point: the matter- 
dominated solution associated with decelerated expan- 
sion, for its DBI-dual j{X) there are found: i) the infia- 
tionary de Sitter late-time attractor (critical point Pds), 
ii) the - also inflationary - de Sitter twin state (point 
Ptw)', a saddle in the phase space, and, in) the point P* 
- associated with a decelerated expansion-phase - where 
the non-linear curvature effects mimic radiation. The 
latter being also a saddle point in ^t. 

Although for other values of the parameter n, the ef- 
fects of the DEI modification are not so spectacular as 
for n = 1, nonetheless, these are appreciable. In particu- 
lar, i) the co-existence of two inflationary future attrac- 
tors (equilibrium points Pds and Ptw) for given values 
of the free parameter n, and, ii) the chance to explain, 
in a united picture, primordial and late time inflation 
as originated by the non-linear effects of the curvature, 
otherwise. 



VI. DISCUSSION 

It is known since long ago, that quadratic (higher-order 
in general) modifications of general relativity are plagued 
by instabilities. Amongst them, we can name, the Ostro- 
gradski, Ricci and tachyon instabilities, and the presence 
of Weyl ghosts - also known as poltergeist, etc. To the list 
here we added the requirement of non-negativity of the 
effective gravitational coupling upon linearization, a sub- 
ject that has not been much discussed in the literature, 
but, as we have shown in subsection A of section II, is of 
not less importance than the other kinds of catastrophic 
diseases of quadratic theories. Actually, linearization of 
an /(AT, ....)-theory around vacuum, maximally symmet- 
ric spaces of constant curvature Xq, leads to multiplica- 
tion of the coupling constant 1/2k^ by an overall factor 
ao = — XqJ'^^, which, in general, is not restricted to 
be positive. For the particular case of an f{X) theory, 
Ricci stability, absence of tachyon instability, and non- 
negativity of the effective gravitational coupling, are not 
independent requirements. Hence, for instance, a Ricci 
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stable /(X)-theory with positive effective gravitational 
coupling upon linearization, is also free of tachyon insta- 
bility. Usually, for most f{X, ...)//(X)-tlieories found in 
the literature, not all of the above mentioned instabil- 
ities are surmounted at once. For instance, f{X,Y,Z) 
theories where the invariants Y and Z enter in the com- 
bination Z — 4Y (like in the Gauss-Bonnet invariant), 
are poltergeist-free, while their linearization is also Os- 
trogradski stable. However, in general, these are either 
Ricci unstable, or have a scalar tachyon instability, or 
both. This is not to speak about positivity of the effec- 
tive gravitational coupling upon linearization. Strictly 
speaking, most one can expect is to avoid really catas- 
trophic instabilities such as the Ostrogradski and Ricci 
ones, and, after our results, also to allow for a positive 
definite effective gravitational coupling. Hence the ques- 
tion: does further modification of quadratic theories re- 
ally matter?, in particular the one that gives name to 
the present paper: does Dirac-Born-Infeld modification 
of quadratic theories really matter? 

The results discussed in this paper suggest that, al- 
though the DBI deformation does not affect the Ostro- 
gradski stability, other important instabilities such as the 
Ricci and tachyon ones are indeed surmounted after ap- 
plying the DBI modification procedure to the original 
theory. Sometimes this is achieved at the cost of re- 
nouncing to the original motivation of the DBI strat- 
egy itself: to avoid singularities. Next one has to care 
about a by-product of the DBI deformation: its impact 
on the cosmic dynamics. It is for sure that the DBI mod- 
ification also affects the dynamics, since the replacement 
/ — > eA(-\/l + 2e//A— 1), inevitably affects the field equa- 
tions. In terms of variables of the phase space this is quite 
clear: as a first visible effect, due to the introduction of 
a new energy scale A, the DBI procedure increases the 
dimension of the phase space. ^ Hence, the question now 
is: would the effect of the DBI modification broaden the 
possibilities of a given model to do cosmology, or would 
it favor the contrary effect? To answer to this question, 
since finding exact solutions to the modified cosmological 
equations is bizarre difficult, one can rely on the tools of 
the dynamical systems. In this regard, knowledge of the 
equilibrium points in the phase space corresponding to a 
given cosmological model is a very important information 
since, independent on the initial conditions chosen, the 
orbits of the corresponding autonomous system of ODE 
will always evolve for some time in the neighborhood of 
these points. Besides, if the point were a stable attrac- 
tor, independent of the initial conditions, the orbits will 



Although, for the example studied in this paper this is the case, 
in general, this is not always true. For instance, for the model 
studied in [sSl the dimension of the phase space is unchanged 
under the DBI replacement. However, even if the dimension is 
not affected, the topology of the phase space is in fact modified 
through, for instance, the occurrence of bifurcations in the space 
of parameters. 



always be attracted towards it (either into the past or 
into the future). Going back to the original cosmologi- 
cal model, the existence of the equilibrium points can be 
correlated with generic cosmological solutions that might 
really decide the fate and/or the origin of the cosmic 
evolution. In a sense the knowledge of the asymptotic 
properties of a given cosmological model is more relevant 
than the knowledge of a particular analytic solution of 
the corresponding cosmological equations. While in the 
later case one might evolve the model from given initial 
data giving a concrete picture that can be tested against 
existing observational data, the knowledge of the asymp- 
totic properties of the model gives the power to realize 
which will be the generic behavior of the model without 
solving the cosmological equations. In the dynamical sys- 
tems language, for instance, a given particular solution 
of the Einstein's equations is just a single point in the 
phase space. Hence, phase space orbits show the way the 
model drives the cosmological evolution from one partic- 
ular solution into another one. Equilibrium points in the 
phase space will correspond to solutions of the cosmolog- 
ical (Einstein's) equations that, in a sense, are preferred 
by the model, i. e., are generic. The lack of equilibrium 
points that could be correlated with a given analytic so- 
lution of the model, amounts to say that this solution is 
not generic, otherwise it can be attained under a very 
carefully arrangement of the initial conditions only. 

As an example, to illustrate the spectacular effect of 
the DBI deformation procedure on the cosmic dynamics 
generated by arbitrary gravity theories, let us consider 
general relativity, i. e., we take the Lagrangian density 
2K^Cg = yJ\g\X, which is a particular case of the theory 
f{X) = ^X"- studied in subsection B of section V. Before 
going into the details of the modification of the asymp- 
totic properties of this theory, by the effects of its DBI 
deformation, we have to briefly discuss about the sta- 
bility of the linearization of the theory, around vacuum, 
maximally symmetric spaces of constant curvature, since 
this topic is central both, to discuss about the physical 
content of the theory, and to study its weak field limit. To 
start with, while the massless graviton is the only propa- 
gating linearized degree of freedom in general relativity, 
the linearized gravitational spectrum of its DBI-dual 



consists of the massless graviton plus a spin-0 massive 
propagating mode. We know, however, that this is also 
an achievement of f{X) modifications of general relativ- 
ity, so that, this was, in fact, an expected effect. To 
see the advantage of DBI-/(X) = J{X) over just f{X), 
consider the simplest "power-law" modification of gen- 
eral relativity f{X) — /iX" (/i,n > 0). When linearized 
around vacuum, maximally symmetric spaces of constant 
curvature Xq, the mass squared of the spin-0 propagating 
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mode is given by 



or 



Xo 
3 



1 



while the Ricci stability requirement, and positivity of 
the effective gravitational coupling, amount to 



and 



- Xo/^' = -n(n - 2)mXo"-i > , 



/^' = n(n-l)MXo"- 



respectively. One sees that: i) for < n < 1 the model 
is Ricci unstable, shows scalar tachyon instability, and 
is unphysical, due to negative effective gravitational cou- 
pling, ii) for 1 < n < 2 the model is stable and physical, 
while iii) for n > 2, it is Ricci stable, but develops scalar 
tachyon instability, and is unphysical due to negative cou- 
pling. Hence, the only physically meaningful range of the 
free parameter n is 1 < n < 2. For its dual DBl- f {X) , 
or f{X) theory, one can check that the mass squared of 
the scalar degree of freedom upon linearization can be 
written as follows: 



mn 



Xo 
3 



(n - 1)A + e(n - 2)fiX'„' 



while Ricci stability and positivity of the effective gravi- 
tational coupling, are now given by 



(n - 1)A + e{n - 2)/xX^' > , 



and 



(n - 2)A + e(n - 4)/iX^ < , 

respectively, with the simultaneous requirement that 

2e^Xj -t- A > . 

For the e = -|-1 DBI embedding, the latter bound is al- 
ways satisfied, while for e = — 1, it is satisfied whenever 
X^ < A/2^. For e = -1-1, the /(X)-theory is Ricci sta- 
ble and physically meaningful (positive definite effective 
gravitational coupling) if, either 



or 



xs< 



xs> 



l\ A 



2 — n J jj, 



1 < n < 2 



4 - n 



-, 2 < n < 4 . 



Meanwhile, for e = —1, it is both, Ricci stable and phys- 
ically meaningful if, either 



1 



n 



-<xs< 



-, < n < 1 



X^ < {^) -, 1< n< 2 
\4-nJ n 



We see that, in both cases: e = ±1, the range of the pa- 
rameter n where the the theory is Ricci stable and phys- 
ically meaningful (and consequently, also free of tachyon 
instability), is wider in the DBI-dual ,f(X) with respect 
to the f{X) modification of general relativity. 

Now we are in position to discuss possible modifica- 
tions of the cosmic dynamics by the DBI procedure. Go- 
ing back to general relativity f{X) = X, the only equi- 
librium point in the phase space, is the one associated 
with matter dominance, where the peace of the cosmic 
expansion is decelerating. The phase space structure of 
its DBI-dual J{X) = e\{.JT+2^XjX - 1), is in fact, 
more complex and rich. Actually, according to the re- 
sults presented in subsection C of section V, there are 
found 3 equilibrium points: i) the inflationary de Sit- 
ter attractor P^S = (2,-1,-4) (0^ = 0, Wcur = — 1, 
q = —1), ii) the "twin" inflationary (de Sitter) solution 
Ptw = (2,0, -(5 - \/5)/2), and iii) the "radiation" -like 
critical point Prad = (0,1,0) (17„, Wcur = 1/3, q = 1), 
associated with decelerated expansion of the universe. 
The last two points Ptw and Prad, are saddle equilibrium 
points in the phase space. An interesting feature of the 
modifled / (X)-theory, is associated with the existence of 
the inflationary (de Sitter) solutions Pas, ^nd Ptw Both 
differ in their Hubble rates of expansion. Actually, since 
fjbz = eA(6iI^)~", then, for Pas we have that 



while, for Pt, 



eX 
24^ 



eA 



3(5 - v^)^ 



Notice that the above solutions exist only for the e = — 1 
DBI embedding. It is clear that, thanks to the existence 
of Pds and Ptw, there is room in the modified ,f{X)- 
theory to explain, in a united way, primordial inflation 
(saddle equilibrium point Ptw), and the present stage of 
the cosmic speedup (late-time attractor Pds)- According 
to this picture, both inflationary stages have their ori- 
gin in the non-linear curvature effects. The fact that the 
primordial inflation is associated, in this case, with a sad- 
dle equilibrium point, explains in a natural way the exit 
from inflation as due to the (in)stability properties of the 
corresponding solution. We have to underline, notwith- 
standing, that this is just a toy model that might not be 
capable of accommodating the existing amount of obser- 
vations. In particular, this objection is obvious in what 
concerns to the formation of structure, since there are no 
equilibrium points that could be associated with matter 
dominance. 
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VII. CONCLUSIONS 

In this paper we have performed a thorough investi- 
gation of the effects of applying the Dirac-Born-Infeld 
strategy to quadratic modifications of gravity. We paid 
special attention to the subtleties associated with lin- 
earization of these theories, as well as of their DBI-duals, 
around vacuum, maximally symmetric spaces of constant 
curvature, an issue that is central to discuss about the 
weak field limit, as well as on the physical content of these 
theories. We have shown, in particular, that one aspect 
to take into consideration under the above linearization 
procedure, is the check of the positivity (non-negativity) 
of the effective gravitational coupling. A theory whose 
linearization fails to give a positive (effective) gravita- 
tional coupling, is devoid of physical meaning. For the 
particular case of f{R) theories, Ricci stability, absence 
of tachyon instability, and positivity of the effective gravi- 
tational coupling upon linearization, arc not independent 
requirements. It has been demonstrated, also, that the 
severe instabilities inherent in quadratic modifications of 
general relativity, can be surmounted - at least smoothed 
out - by considering their DBI-duals. This achievement 
is, in some cases, at the cost to renounce to the original 
motivation of the DBI deformation strategy: to remove 
singularities. 

An important focus of this paper has been the study 
of the modifications carried on the asymptotic properties 
of quadratic gravity by the DBI procedure. By applying 
very simple recipes of the dynamical systems, it has been 
demonstrated that the structure of the phase space is, in 
fact, modified by the DBI deformation strategy. The di- 
mension of the phase space is increased, which means, 
in turn, a richer asymptotic structure. Even if, in many 
cases, the results are not spectacular, we think that con- 
sidering DBI-duals of quadratic modifications of gravity 
can be an interesting arena where to look for alterna- 
tive explanations of such important cosmic mysteries as 
primordial inflation and the present cosmic speedup. 
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VIII. APPENDIX: DYNAMICAL SYSTEMS 

Here we include brief tips of how to apply the dynami- 
cal systems tools in general. In order to apply these tools 
one has to follow the steps enumerated below: 

1. To identify the phase space variables that allow 
writing the system of cosmological equations in the 



form of an autonomus system of ordinary differen- 
tial equations (ODE), say:^*^ 

2. With the help of the chosen phase space variables, 
to build an autonomous system of ODE out of the 
original system of cosmological equations (r is the 
time-ordering variable, not necessarily the cosmic 
time): 

dxi 

= fi{Xl,X2,-Xn) . 

Notice that the RHS of these equations do not de- 
pend explicitly on r (that is the reason why the 
system is called autonomous). 

3. To identify the phase space spanned by the chosen 
variables (xi,X2, ...a;„), that is relevant to the cos- 
mological model under study. This amounts, basi- 
cally, to define the range of the phase space vari- 
ables that is appropriate to the problem at hand: 

* = {{xi,X2, ---Xn) : bounds on the Xi-s} . 

4. Finding the equilibrium points of the autonomous 
system of ODE, amounts to solve the following sys- 
tem of algebraic equations on {xi,X2, ■■■Xn)- 

Mxi,X2,:.Xn) = . 

5. Next one linearly expands the equations of the au- 
tonomous system of ODE in the neighborhood of 
the equilibrium points pk = Pk{xi,X2, ■■■Xn), k = 
1,2, ...m'}^ I. e., one replaces Xi ^ Xi + e^, where 
Ci are the small (linear) perturbations around the 
equilibrium points. Hence the system of ODE be- 
comes a system of linear equations to determine the 
evolution of the Cj-s: 

otherwise, since /j = fi{p) = 0, then 

3 

where we have introduced the linearization or Ja- 
cobian matrix [Ml] = dfi/dxj. 



There can be several different possible choices, however, not all 
of them allow for the minimum possible dimensionality of the 

phase space. 

In general the number of equilibrium points is different from the 
dimension of the phase space; m ^ n. 
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6. The next step is to solve the secular equation to de- 
termine the eigenvalues of the linearization ma- 
trix at the given equilibrium point p: 

det|M(p)^-A[//| = 0, 

where [J7/] is the unit matrix. 

7. Once the eigenvalues of the linearization around a 
given equilibrium point p have been computed, the 
evolution of the perturbations is given by 

n 

e^ir) = ^(eo)^exp(AjT) , 

j 

where the amplitudes (cq)^ are constants of inte- 
gration. 

If all of the eigenvalues have negative real parts, the per- 
turbations decay with r, i. e., the equilibrium point is 



stable against linear perturabtions. The corresponding 
equilibrium point is said to be a future attractor. If at 
least one of the eigenvalues has positive real part, the 
perturbations grow with r so that these are not stable in 
the direction spanned by the given eigenvalue. Hence the 
point is said to be a saddle. The perturbations around 
a given equilibrium point are unstable, in other words 
the point is a past attractor (a source point in the phase 
space), if all of the eigenvalues have positive real parts. 
Points whose linearization is characterized by complex 
eigenvalues are said to be spiral equilibrium points, and 
are commonly associated with oscillatory behavior of the 
corresponding solution. If at least one of the eigenval- 
ues has a vanishing real part, the equilibrium point is 
said to be non- hyperbolic. In the latter case, in general, 
and unless some of the non-vanishing real parts of the 
eigenvalues are of opposite sign, one can not give conclu- 
sive arguments on the stability of the equilibrium point. 
Other techniques have to be applied. 
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